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Gaussianity and statistical isotropy of the universe are modern cosmology's minimal set of hy- 
potheses. In this work we introduce a new statistical test to detect observational deviations from 
this minimal set. By defining the temperature correlation function over the whole celestial sphere, 
we are able to independently quantify both angular and planar dependence (modulations) of the 
CMB temperature power spectrum over different slices of this sphere. Given that planar dependence 
leads to further modulations of the usual angular power spectrum Ci, this test can potentially reveal 
richer structures in the morphology of the primordial temperature field. We have also constructed an 
unbiased estimator for this angular-planar power spectrum which naturally generalizes the estimator 
for the usual Ce's. With the help of a chi-square analysis, we have used this estimator to search for 
observational deviations of statistical isotropy in WMAP's 5 year release data set (ILC5), where we 
found only slight anomalies on the angular scales £ = 7 and ^ = 8. Since this angular-planar statistic 
is model-independent, it is ideal to employ in searches of statistical anisotropy (e.g., contaminations 
from the galactic plane) and to characterize non-gaussianities. 

PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

Many efforts have been made towards understanding the statistical properties of the cosmic microwave background 
(CMB) temperature field in the past few years. The main motivation behind these efforts is that, in a homogeneous 
and isotropic universe in which inflation is driven by a single canonical scalar field, the primordial temperature field 
is set by Gaussian and statistically isotropic physical processes. Since nonlinear evolution destroys all these putative 
initial gaussianities, we must search for any fundamental deviations from these statistical properties at early epochs 
and as close as possible to the linear regime. This makes the CMB the ideal physical observable to employ in searches 
of statistical anisotropies and non-gaussianities. Any significant observational deviation from this picture could reveal 
us something as yet unsuspected about the basic nature of our universe. 

While this program seems to be well motivated by itself, careful analysis of recent temperature maps obtained by 
the WMAP team |T1|2]|3] have hinted at some apparent anomalies - mainly in the low multipoles sector jH |5| |6j [TJ |8] . 
K one leaves aside for a moment the perennial problem of a posteriori statistics j^, these findings raise the possibility 
that the anomalies could be a first hint towards some new physics. Theoretical attempts to explain their origin 
include primordial magnetic fields [10], non-trivial cosmic topologies jlll 112) . globally anisotropic models of the early 
universe |131 [Til [TSl [T6| as well as local manifestations of cosmic anisotropy [TJ, and even anisotropic models of 
dark-energy |181 [TO] . Of course, there is a good chance that these anomalies are due to astrophysical effects [20] or 
even some residual instrumental cross-contamination, in which case our universe can still be easily accommodated in 
the standard scenario. It is a question of utmost concern to decide whether these known anomalies (as well as others 
which may be found in the future) are isolated statistical flukes, or if they are due to new physical/astrophysical 
effects. 

Despite its importance and the efforts spent on it, we still have no compelling explanation for the nature of the 
\ow-£ anomalies. The main difficulty is twofold: first, we still do not know how to optimally separate the question of 
gaussianity from that of statistical isotropy (see however [21] for a first step in this direction.) It is thus possible that 
our universe is Gaussian but statistically anisotropic, statistically isotropic but non-gaussian or even non-gaussian and 
anisotropic. Second, if the universe is neither Gaussian nor statistically isotropic, then it can be - from a statistical 
point of view - virtually anything: there is only one kind of gaussianity and isotropy, but there are infinite ways to 
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brake either one. The absence of theoretical guidelines will inevitably lead to an infinite number of models and no 
underlying symmetries, which still would mean that we could not account for confirmed anomalies. 

This means we must analyze the problem in as much a model-independent manner as possible. We can, for example, 
start from the very basic definition of our statistical quantities (such as the two-point correlation function) and check 
whether they can be modified in a model-independent manner, basing our reasoning solely on physical symmetries 
and observational hints. 

One such possibility is to consider the two-point temperature correlation function, C(ni,n2), without some of the 
symmetries of the underlying space-time. Attempts in this direction have been made by Pullen and Kamionkowski 
|22) . where the temperature correlation function is assumed to depend on the direction of any given unit vector in 
the celestial sphere, in such a way that one can search for power multipole moments in temperature maps. However, 
that approach consists of considering the temperature correlation function at zero lag, and thus it does not allow us 
to consider correlations between two different points in the sky. 

Another possibility is to consider the correlation function in its full form, i.e, a function that depends on all pairs 
of independent unit vectors in the sphere 5^. This idea, which was introduced by Hajian and Souradeep |23l l24l [25] . 
consists in expanding the temperature correlation function in a bipolar spherical harmonic series in order to take into 
account its functional dependence. The authors then construct a bipolar power spectrum ki which can account for 
deviations from statistical anisotropy if observations give us K£ > at a statistically significant level. Unfortunately, 
that approach is too generic: it is not clear what the associated statistical test is measuring, or how one can motivate 
it in the absence of an underlying theoretical or phenomenological model. 

In this work we also go back to the two-point correlation function C(ni, 712), but instead ask whether it can depend 
not only on the separation angle between two given unit vectors, cosi? — hi ■ 1x2, but also on the orientation of the 
plane of the great circle defined by the unit vectors. Such a functional dependence can be unambiguously constructed 
once we realize that for any two unit vectors in the CMB sky, their angular separation and their associated plane are 
uniquely defined by their dot and cross products, respectively. The new planar dependence (on the direction defined 
by the normal to the planes of the two unit vectors) codifies modulations of the usual two-point correlation function 
as we rotate these planes while keeping the separation angles ■& fixed. 

We have also constructed, in a completely model-independent way, an angular-planar power spectrum and its 
associated unbiased estimator, which naturally generalizes the usual angular power spectrum Ci, and for which we 
recover the known results in the limit of statistical isotropy. 

Our approach has a strong observational motivation, which lies in the fact that some astrophysical planes, 
like the galactic and ecliptic ones, play an important role in CMB measurements and could still be manifested 
in the data if the foregrounds were improperly removed. One such example was possibly found in [26] where, 
besides the alignment of the multipoles £ = 2 and £ = 3, the authors detected a strong correlation between 
these two and the ecliptic plane. The existence of a preferential plane could also be related to the so-called 
north-south asymmetry [5] |8] [53 , in which case a plane could naturally separate regions of maximum and minimum 
temperature power. There exists also a third situation in which a physical plane can play and important role in 
cosmology, namely, the unavoidable presence of our galactic plane in all CMB measurements acts as an important 
source of astrophysical and foreground contamination. All these facts lead us to believe that a planar signature 
on the correlation function would be an important statistical property of the CMB, and is a potential test of its nature. 

We have organized this work in the following way: we begin ^jlljwith a brief description of the two-point correlation 
function and its general properties. After discussing some of its known generalizations, we extend our argument to 



include a planar dependence. In ^III we carry a multipolar decomposition of the correlation function with planar 
dependence and show how the resulting coefficients (i.e., the angular-planar power spectrum) are related to the usual 
temperature multipolar coefficients a^m's. This leads us to the question of how to build an unbiased estimator to 
measure planar signatures in temperature maps and, in particular, how this can be implemented with the help of 
a simple chi-square analysis. We illustrate, still in this section, the application of our statistics to the well-known 
ACDM concordance model, where we present some figures for the "best-fit" ACDM angular-planar power spectrum. 
In |IV| we use a chi-square test to search for planar signatures in the WMAP full-sky temperature maps, and show 
that the angular scales i — 1 and ^ = 8 seem to be slightly anomalous for a particular range of planar separation 
We conclude in ^Vl where we also give some perspective of further developments. 



II. TEMPERATURE CORRELATION FUNCTION 



The main observable in the CMB is the temperature fiuctuation field, AT. In its full generality, this field is a 
function of a position vector n and of the time interval in which we measure this temperature - but in practice our 
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measurements are made in time intervals which are negligible compared with the cosmological timescales. The field 
AT is a scalar, continuous function on the unit sphere, which means we can decompose it in the usual fashion, in 
terms of spherical harmonics: 

AT(n) = ^a,„,r,„(n). (1) 

i,m 

All information is therefore encrypted in the multipolar coefficients aim- Essentially all infiationary models predict 
these coefficients not as uniquely given, but rather as realizations of a random variable, in such a way that the 
physics is not in the aim's themselves, but rather on their statistical properties. Since, by construction, this field has 
zero expectation value, (AT) — 0, the two-point correlation function expresses the first nontrivial momenta of the 
underlying statistical properties of the physical field, and is given by: 

C(ni,n2) ^ (AT(ni)AT(n2)> = E E («^i™,<™.)y^.m.("i)>^£:™.(^2) • (2) 

£i,mi £2, "12 

Alternatively, the covariance matrix above, {o-iimiO-i^^^) , gives all the information about the quadratic momenta of 
the underlying distribution. If the field AT is Gaussian, then this covariance matrix encloses all the information 
that is needed to describe the nature of the fiuctuation field ([ij. In this work we shall restrict ourselves to a fiducial 
Gaussian model, for simplicity. 

We note also that the separable nature of the definition ([2} implies a reciprocity relation for the correlation function: 

C(ni,n2) = C(n2,ni). (3) 
This symmetry must always be satisfied, regardless of the underlying physics. 



A. Isotropic case 

In a globally homogeneous and isotropic universe, the two-point correlation function of the temperature can only 
depend on the separation angle between the vectors hi and n2, that is: 

C{ni,fi2) = C{hi ■ h2) 

--^CiPeihi-h^). (4) 
47r 



E 



Comparing this expression with Eq. ([2}, we notice that the covariance matrix becomes diagonal: 

{O'iimial^m^) — Gi^bl^l^^mxrai , (5) 

with the diagonal terms given by the angular power spectrum, Ci. In principle the angular power spectrum suffices to 
describe the statistical properties of the temperature field ([T|. However, since we have only one universe to measure, 
and therefore only one set of aim's, the average in ^ is poorly determined. The best we can do then is to take 
advantage of the ergodic hypothesis, which states that averaging over an ensemble can be treated as averaging over 
space, and hence to consider each of the 2£ + 1 real numbers in agm as statistically independent, in such a way as to 
build a statistical estimator for the C^'s: 

1 ^ 

m—~£ 

Since (C/) — Ci, this estimator is said to be unbiased. Also, because for a Gaussian field {{C i — C i){C ii' — C i')) oc (5«', 
this estimator has the least "cosmic variance". Ci is, therefore, the best estimator that can measure the statistical 
properties of the multipolar coefficients aim when both statistical isotropy and gaussianity hold. 



B. Some anisotropic cases 

The first line in Eq. Q for the temperature two-point correlation function is valid if and only if the universe 
is statistically isotropic. This means that any functional dependence that does not reduce to a dependence on 
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cos ?9 = ni • Ti2 will measure some deviation from statistical isotropy. There are infinite possible combinations of fii 
and TJ,2 that violate statistical isotropy. However, since the vectors hi and n2 are constrained to have a common 
origin and size, symmetry and simplicity does not leave us many choices. One possibility is to consider these two 
vectors as being the same, in which case we are left with a correlation function of the form: 

C : 5*2 ^ M , (6) 

and for which a decomposition similar to exists. This form of the correlation function makes it suitable for searching 
for power multipole moments in CMB temperature and polarization maps, once we define a power multipole moment 
estimator [22]. On the other hand, this is also a correlation function at zero lag, so by construction it does not allow 
us to consider anisotropic correlations between different points in the sky. 

A second possibility is to consider the correlation function as being the most general (but separable) function of 
two unit vectors that one can possibly have |23) : 

C-.S^xS'^^R. (7) 

This function admits a decomposition in terms of the bipolar spherical harmonics |23] which has the nice property of 
behaving - in many mathematical aspects - as the usual spherical harmonics. The main drawback of the decomposition 
([7]l, however, is that it carries too many degrees of freedom which, in the absence of a specific cosmological model, 
cannot be resolved with simple estimators. Therefore, these two approaches are either too simple or too generic to 
reveal deviations from statistical isotropy in a more model-independent way. 



C. Anisotropy through planar dependence 

The guiding principle used in the construction of ([gJ and Q is rather general and based mainly on our prejudices 
about what statistical anisotropy should look like. However, in the absence of theoretical guidelines, we have to confine 
ourselves to the observations of the CMB temperature or, more specifically, to the signature of its known anomalies. 
One example is the role played by the galactic and ecliptic plane in the quadrupole-octupole/north-south anomalies 
m H] m [71 18], not to mention the importance of our galactic plane as a source of foreground contamination in the 
construction of cleaned CMB maps. The existence of a cosmic plane might even be a manifestation of some mirror 
symmetry |28) . 

In general, the simple fact that we are bound to make all our measurements inside our galactic plane suggests that 
the correlation between fields at two positions fii and n2 might be sensitive not only to their separation angle but 
also to the orientation of the plane they live in, as is shown in Fig. [l] 




Figure 1: Schematic representation of the functional dependence of the correlation function (|8|. In the anisotropic case we are 
considering, fields at the positions given by vectors living on different planes can have different correlations, regardless of their 
angular separation. 

Such a planar dependence can be included in the correlation function if we realize that two unit vectors on the 
sphere S'^ uniquely define both a separation angle and a direction n perpendicular to the great circle (or plane) 
where they live. We are then left with a new possibility for the functional dependence of the two-point correlation 
function: 



C(ni,n2) = C{ni x 72.2) , 



(8) 
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which corresponds formally to a function of the form C : M, where is the set of all points {x, y, z) such that 

+ 1p + < 1 IIO]. Defining n = ri.i x n2, the above expression can be further decomposed in spherical coordinates 
as follows: 

C(n) = VV?^C^P,(cosz?)yz™(n), I e 2N , (9) 

where: 

|n| = sini? , n = {9, (j)} . 

Notice that: 

hi ■ h2 = cos 7? = cos 9i cos 02 + sin 6i sin 62 cos(0i — (j>2) , (10) 

where {0i,(j)i) are the angles defined by the vectors n^. 

Some comments on the decomposition ([9} are in order. First, we note that there is an intrinsic ambiguity in the 
sense of the vector n (as we might as well have defined n = n2 x fii), which is obviously inherited from the ambiguity 
in the definition of the normal to a plane. This ambiguity can be avoided if we restrict the sum in / to even values, 
which is what we will do from now on. Note that such a restriction arises naturally as a consequence of the reciprocity 
relation Eq. ([3}. Second, for £ = we recover ([gJ and therefore all the analysis made in j22] arises as a special case 
here. 



III. ANGULAR-PLANAR POWER SPECTRUM 

The multipolar C^"' coefficients in Eq. ^ correspond to a generalization of the usual angular power spectrum 
C^'s. In fact, they can be seen as a spherical harmonic decomposition of the angular power spectrum, if it suffers 
modulations as we sweep planes on the sphere. The function Ce{n) for a given £ is: 



C,(n) = ^/4^^C^'i^„(n), I e 2N . 



-'00 



Clearly, the monopole of Ce{n) (the average over the whole sphere) is the usual angular power spectrum, — Ci, 
and the higher multipoles measure modulations of the spectrum. 

Since we are restricting our analysis to the Gaussian case, the set of coefficients C[™ completely characterizes the 
two-point correlation function. Still, what is accessible through observations are temperature maps which we can use 
to try to estimate the correlation function. In this respect the multipolar coefficients C^™ would be of limited interest, 
unless we can relate them directly to our observables. It would be interesting if we could, for example, relate these 
coefficients to the covariance matrix {cL£-^miO,e^m2^ equating expressions (|9jl and ([2|), as is usually done. However, 



this procedure is far from being trivial, since the complicated coupling of the angles i?, 9 and <j) defined in (10 1 make 
it difficult to use the usual orthogonality relations to isolate the C^""s. 

Fortunately, as we show below, we can estimate the C]™'s if we use the invariance of the scalar product ni • h2 and 
chose our coordinate system in order to integrate out the d dependence. Once this is done, we make a passive rotation 
of the coordinate system and then we integrate over the remaining angles 9 and 0, which then are given precisely by 
the Euler angles used in the rotation. The details are rather technical and can be found in the Appendix. The final 
expression is: 



where the 6-index expression in parenthesis is the Wigner 3J symbol, and 

^E(-ir^^i'"^^^'" ( J t -m) ^<-cos^)P,(cos^)e^™^ (12) 

where Xi-m are a set of coefficients resulting from the 1} integration, which vanish unless ii + m — even (see the 
Appendix for more details.) 
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It is easy to show that expression ( [TT| induces no coupHng between the eigenvalues £ and I, as expected, since the 
length of the vector n is completely independent of its orientation. There are, however, subtle couplings present in 
11 1 which do make a difference when we apply it to real data. This is due to the Legendre polynomial in the integral 
12 I, which selects only those values of £i and £2 which have the same parity as the angular momentum £. Moreover, 



the 3 J symbols appearing in (111 give different weights to the triple {I, £1,12) depending on the parity of {£i,£2) 



and, as a consequence, we can expect typical oscillations in any function of (111 that we may build when plotted as 



a function of £. This will be shown explicitly in the next section, when we apply these tools to the WMAP 5 year data. 



Expression (111 does not take into account the fact that real data is not given exactly by ([ij, but rather by 
a pixelized temperature map which is a combination of the true cosmological signal, plus instrumental noise and 
residual foreground contamination. Schematically, the temperature of the map in each pixel i is given by 

Typically, the cosmological signal AT"^ is smoothed out by a Gaussian beam W{fii) of finite width which, in harmonic 
space, is given by Wi = exp(— f^cr^), where ab = 6'fwiim/'\/81n2 and 6'fwhm is the beam full width at half maximum. 
For the V-band frequency map of the WMAP experiment, 6'fwhm — 0.35°, which implies a minimum £ynin ^ 390 for 
which the effect of a beam smoothing will be important, much higher than the low-^ regions where known anomalies 
were reported. Thus, for the sake of simplicity we will neglect the effect of the beam in this work. Also, for the £ < 390 
region, cosmic variance is known to dominate the source of error over instrumental noise, and therefore we can neglect 
the latter as well. On the other hand, the residual foreground can be an important source of contamination, and 
therefore deserves a careful analysis which is beyond the scope of the present work. In a companion paper we carry a 
more rigorous analysis of planar signature in CMB data in which the effect of the residual foreground will be estimated 
I29i. 



A. Statistical estimators and x analysis 



We now would like to use expression (11 1 to examine the observed universe. We start by noting that in the limit 
of statistical isotropy (SI), that is, when (a^jm^a^jmj) = C£-^Si-^i2^rnim2^ expression (111 reduces to 



C 



(13) 



Conversely, if the only non-zero C^™'s are given by I — m — 0, then = Ci. Therefore, statistical isotropy is 
achieved if and only if the C^™'s are of the form ([l3|, and any observational deviation from this relation would be an 



indication of statistical anisotropy. 

However, we only get to observe one universe, and this makes the "cosmic sample variance" a severe restriction that 
we have to live with. This means that if we want to know, let's say, the mean value and variance of the C^""s, we will 
have to build statistical functions which can only estimate these prope rties, just like it is done with the fundamental 
quantities and the associated estimators (see the discussion in ^IIAl. 



In other words, in order to evaluate the statistical properties of the C™'s we will have to treat them as our new 
"fundamental" quantities, which will be determined exclusively as a function of the a^m's. As a consequence, we will 
redefine expression (111 as: 



C 



27rV2Z + 1 



I £1 £2 

m mi 7712 



(14) 



and will treat the coefficients C^™ as uniquely given once we have a map. Of course, expression ( 14 1 is nothing more 



than the unbiased estimator of the angular-planar power spectrum ( |ll| and - as long as cosmic variance is an issue - 
this "second order" approach we are adopting here (i.e, the prescription of adopting this estimator of the correlation 
function as our fundamental quantity, rather than the temperature field) is the best we can do when searching for 
statistical deviations of isotropy. In theory, it is also possible to use the CMB polarization induced by galactic 
clusters to probe different surfaces where CMB photons last scattered, and to use such independent measurements 
as a way to alleviate cosmic variance |30l [31] . However, the gain in terms of a reduced variance is still limited. 

Having these limitations in mind, we can now ask: how good does a theoretical model of anisotropy, Cf^''™", fit 
the observational data once it is given by (14 1? To answer this question we can use the well-known chi-square (x^) 
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goodness-of-fit test, which in our case can be written in the following "generalized" form: 



Im 



m— — l 



in which the scales I and i are seem as independent degrees of freedom, and where cr^™ is just the standard deviation of 



the difference (C^™ - Cf Although this expression can be readily applied to any theoretical model of anisotropy, 
in practice it is better to work with its reduced version: 



2\l 



ixt) 



21 



m— — / 



th,/m 1 2 



lm\2 



(15) 



which is just the chi-square function divided by the 21+1 planar degrees of freedom. 



Expression ( 15 1 will be the starting point of our statistical analysis, which we will pursue in detail in { IV Before we 
move on, it is important to choose a particular cosmological model of anisotropy against which we want to compare 
our estimator (14 1. 



B. ACDM model 



The most important model to be analyzed using the estimator (14 1 is, of course, the concordance ACDM model 



which was confirmed with striking accuracy by the 5 year release dataset of the WMAP team [2]. For this model, 
statistical isotropy holds and, as we have shown, any multipolar coefficient Cl™ with non-zero planar dependence 



should be identically zero in this case [see expression (13 1.] Therefore, for this model we can take: 



(16) 



where it should be clear that we are only considering the cases with / > 2 (as we will do from now on). Now, if the 
data under analysis is really Gaussian and SI, then its covariance matrix can be explicitly calculated: 



J^U'nim' 



iic. 



Ini 



Cf ''™)*(C, 



I'm 



c 



th,l' 



)> 



Sn'S„ 



where we have used ^ and the null hypothesis (16 1. 

This covariance matrix has some interesting properties: first, we note that the planar degrees of freedom in (15 1 
are really independent in this case. Moreover, the variance (cr^™)^ — Al"™™ 



= Stt 



Cf, C'f. 



becomes m- independent: 

2 



(17) 



Second, its diagonal terms (i.e., a\) are completely determined by the angular power spectrum Ci , up to some 
geometrical coefficients which arise as a consequence of the way in which we split our CMB sky. This makes it 
possible to give a visual interpretation of the angular-planar power spectrum C^™, similar to that of the C^'s. For 
that, let us introduce the reduced angular-planar spectrum: 



2i+l 



(18) 



which has a simple interpretation when compared to the usual angular spectrum, because 
as can be easily shown using Eqs. (17 1 and (25 1. 

In Fig. [2] we show some plots of the reduced spectrum Hj, both as a function of / and £. Notice that, as a result 
of our planar splitting of the CMB sky, the low-^ sector of the spectrum Hg is suppressed when we consider planes 
separated by smaller angles (bigger values of I). This is a consequence of the nontrivial coupling of the moments 
£1 and £2- since the C^'s are roughly given by a monotonically decreasing sequence, and since \l — £i\ < £2 ^ I + £1, 
bigger values of I make the moment £2 probe deeper and deeper regions of the Sachs- Wolfe plateau. This suppression 
reaches cosmological scales up to the first acoustic peak, after which the planar dependence becomes negligible. 
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Figure 2: Reduced angular-planar power spectrum for a fiducial ACDM model. The left panel shows the spectrum as 
a function of £ for the particular values (from top to bottom) I = (0,2,4,6,12); the case 1 = (dotted line) represents the 
standard angular power spectrum Ce. The right panel shows the same spectra, but now as a function of I, for the particular 
values £ = (2,3,4,6,8). 



IV. 

We now come back to the question of how the angular-planar power spectrum C^™ fits the observed universe. We 
begin by showing that if we want to compare our data against the standard ACDM universe, then the chi-square 



function (15 I becomes a very simple expression. As we have shown in the preceding section, for this model, a ™ = a\ 



and C*''''"' = 0. Therefore (|15| simplifies to: 



m— — l ^ 



where C? ^' ^ is calculated by applying the estimator (14) to the data given by a^^. It is now clear that if the data 



under analysis is really Gaussian and statistically isotropic, then it should be true that: 

1 ^ / /^obs,/m w^obs Jm\ 

m— — l ^ 

This means that a positive test of planarity will be quantified by how far our chi-square function deviates from unity. 
We can do even better and define a new function as: 

Xl ^ {xl)\ - 1 (20) 

which, if significantly different from zero, will point towards anisotropy. 

It should be stressed that, for a given CMB map, the chi-square analysis must be done entirely in terms of that 



map's data. Indeed, any arbitrary introduction of a fiducial bias in (19 1 (for example, by calculating cr^ using Cf''^^^) 
would only include our a priori prejudices about what the map's anisotropies should look like. The angular spectrum 
Ci, being by construction a measure of statistical isotropy, can only be said to be small/big when compared to a 
particular cosmological model (for example, the ACDM model). Consequently, an anomalous detection of Ci is by 
no means a measure of statistical anisotropy, and it is this value that should be used to calculate a\ if we want to 
find deviations of isotropy, regardless of how high/low it is. Note also that while the function xl has some "isotropy 
variance" which could be computed for the ACDM model from first principles, in practice it is much easier to simulate 
many realizations of a Gaussian and isotropic random field to obtain that variance. 

Finally, we would like to mention that although each number xl is an individual measure of anisotropy (i.e., 
planarity), a consistently biased set of values over a range of Z's or ^'s can also be seen as an indication of anisotropy, 
even if all individual x^'s in that range are well within their variance limits. 



Following the prescription outlined above, we applied the estimator (20 1 to the 5 year WMAP full sky data (also 
known as ILC5 map) where, for practical reasons, we have restricted our analysis to the range of values t g [2, 12] 
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and I G [2, 12] (notice that the momenta I can only assume even values). Our results are presented in Fig. |3] where 
we keep the momenta I fixed, and vary the momenta related to angular separation, £. As discussed before, for this 
range of values cosmic variance dominates over other sources of noise. We estimated the effects of cosmic variance by 
running a simulation of 10'^ realizations of this estimator, using the best-fit (theoretical) scalar Cg's made available in 
|32) : this corresponds to the shaded area in Fig. |3] 

It is also important to explain that, while the data points in Fig. [3] were calculated using the ILC5 map alone, we 
have also included in our analysis a rough estimate of the possible residual foreground contamination present in the 
data. This was done by computing the sample variances of the full-sky maps shown in Table |lj which were then used 
as error bars. In other words, the error bars in Fig. |3] do not account for instrumental noise, which is believed to be 
under control at these angular scales. 



Full sky maps 
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Hinshaw et. al. 
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de Oliveira-Costa et. m 


I. [34] 


Kim et. al. 


[35J 


Park et. al. 


[36J 


Delabrouille et. al. 


[37J 



Table I: Full sky WMAP maps used in our analysis to estimate the error bars of Fig. [3] The data points in this (and subsequent) 
figure was calculated using only the ILC5 map i33j . 



Theses figures present some peculiarities: first, we notice that the magnitude of the error bars oscillate for the 
smallest values of t. As we mentioned in |III| this is partially a consequence of the 3J symbols, which are weights 
appearing in the definition of the anisotropic power spectra and whose effect is to couple differently odd and even 
multipoles. The second peculiarity is that, in all these figures, the modulations of the quadrupolar moment i = 2 are 
entirely consistent with zero. This result suggests that the low value of the quadrupole is perhaps not a consequence 
of statistical anisotropy, at least for the test we are considering here. Note also that the octupole £ = 3, which has 
been reported as unusually planar by some groups, grows slightly from Z = 4 to ^ = 8, although it is compatible with 
cosmic variance in all the planar range considered. 

In what concerns deviations of isotropy, our analysis shows that the most "anomalous" scales are in the sectors 
(Z,^) = (4,7) and (6,8), where we can see that the points x| and x% ^-re only marginally allowed by the 2cr cosmic 
variance area. 
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Figure 3: Anisotropic angular-planar estimator applied to the WMAP ILC5 data. The panel shows Xi as a function of £, for 
the particular values / — (2,4,6,8,10,12). The shaded region and error bars represent 2a cosmic variance and systematical 
(foreground) contaminations, respectively. 



In order to make the visualization of the above figures easier, we repeat the analysis but now keeping the angular 
separation £ fixed and varying the planar separation I. The result is shown in Fig. |4] Notice that the planar 
modulations of the quadrupole £ = 2 are consistently positive, but always compatible with zero. We can also see in 
these figures the growing behavior in the octupole £ = 3 from / = 4 to Z = 8 as mentioned before. 
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Figure 4: Anisotropic angular-planar estimator applied to the WMAP ILC5 data. The panel shows as a function of I, for 
the particular values £ = (2, 3, 4, 5, 6, 7, 8, 9, 10, 11). Note that for fixed I, the error bars do not oscillate (see the text for more 
details). 
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V. CONCLUSIONS 



We have investigated the minimum statistical framework of modern cosmology by enlarging the domain of the 
two-point correlation function to admit not only the usual angular dependence, but a directional (planar) dependence 
as well. Our observable, the anisotropic angular-planar power spectrum, can account not only for the usual angular 
separation between any two spots in a CMB map, but also for any planar signature that this map might have. 
Besides having a strong observational motivation, an interesting feature of this approach is that it leads naturally 
to an unbiased estimator of statistical anisotropy, in the same spirit as is done with the multipolar temperature 
coefficients a^m's and their associated estimator C(. As an example of its use, we have applied this estimator to a 
concrete model of cosmology, i.e., the ACDM model, where we have shown that under the hypothesis of gaussianity 
and statistical isotropy, the angular-planar power spectra have zero mean, but of course, non-zero covariance. 

By means of a simple chi-square analysis, we have also applied our estimator of planar anisotropy to the WMAP 
ILC5 data, where we found that the planar modulations of the quadrupole £ — 2 are compatible with the null 
hypothesis over the range of planar momenta I G [2, 12] we probed. Our results suggests that the low value of the 
quadrupole C2 is perhaps due to some local physics, and not to deviations of statistical isotropy, at least as far as 
planar modulations are concerned. 

Our analysis has also shown that the angular scales 1=7 and £ = 8 suffer some degree of modulation around the 
planar scales I = 4 and I = 6, respectively. This could be an indication of some foreground contamination coming 
from a planar region of typical size A? = 4 ~ 6. However, a complete treatment of the sources of errors and the effect 
of masks is needed before we can reach a more definitive conclusion - for that analysis, see |29] . 

From a theoretical perspective, our techniques can be readily applied to any particular model of inflation predicting 
a specific anisotropic shape for the matter power spectrum. Due to the generality and simplicity of our formulas, the 
angular-planar power spectrum can also be used to analyze CMB polarization. Other possible applications include 
stacked maps of cosmic structure, such as the galaxy cluster catalog 2Mass |38] . 

We finally mention that, although in this work we have focused on testing isotropy while keeping within the Gaussian 
framework, our tools can also be used to search for deviations from gaussianity in a completely model-independent 
way. 
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APPENDIX A 
1. Derivation of ( [iT] ) 

We will present here the details of the derivation of expression ( [TT| . We start by equating expressions ([9| and ^ 

EE^7^^^^^(^°^^)^'"(")== E E («^^™i<™2>^^i™i("i)>^^2™.(^2). (21) 

As mentioned in the main text, the inversion of C^™ as a function of the a^m's is not a trivial task, since the vector n 
depends non-linearly on the vectors hi and 722- The easiest way to achieve this goal is to pick up a coordinate system 
where only the i? dependence (i.e., the modulus of the vector n) is present. After integrating it out, we rotate our 
coordinate system using three Euler angles to recover back the {0, (j)) dependence, which can then be integrated with 
the help of some Wigner matrices identities. We therefore start by positioning the vectors hi and n2 in the xy plane. 



i.e, we chose hi = (7r/2,0i), h2 = {n/2,4>2). By (lOl we then have cosi? = cos((/>i — (f>2)- Using the relation |39] 



y 47r (£+m)!! (£— m)!! 



F,„(V2,</.) = A,,„e™^ where X^m ^ V ' V 4^ (^+")'! (^-™)!! - ^ - ^ ^32) 

otherwise 
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and n = (sin t?, 0,0), we can integrate the i? dependence on both sides of (21 1. This gives us 



'TT 

l.m ii ,mi £2 ,m2 



where we have introduced the following definition 

iM.m, = -\l,mAl.m2 I Pi{cOs{c^l " ^2) )e'('"^ "^^^^^ d(cos(01 - ^2)) • (24) 



We need now to integrate out the 9 and (j) dependence in the right-hand side of (23 1 which was hidden due to 
jr choice of a particul 
jtation of our coordins 
and aim's according to 



our choice of a particular coordinate system. In order to do that, we keep the vectors fii and n2 fixed and make a 



rotation of our coordinate system using three Euler angles lo = {a,/3,7}. This rotation changes the coefficients C]™'s 



aim - E ^mra' {^)o.lm' , C^™ - ^ Dmm' {^)C\"^ 

m' m' 

where C^™ and are the multipolar coefficients in the new coordinate system and where D\^^,{uj) are the elements 
of the Wigner rotation matrix. The advantage of positioning the vectors hi and n2 in the plane xy is that now the 
angles 9 and </) are given precisely by the Euler angles (3 and 7, regardless of the value of a 

Ec^ii™(o,o) = ^ci-' \ y^DLm-{o^,p,i)yUQ,m = Y.c'r'Yi,^m{p,i) 

l,m l.ni' \ m / ;,m' 

where in the last step we have used Yim{0,0) — \/ (2Z + l)/47r (^^o- Therefore, in our new coordinate system we have 
(dropping the "~" in our notation) 

We may now isolate C^™ using the identities |39] 



Stt 



2 



duo D ^* ,(uj)D'^ , (uj) — SKl^Sm.rn^Sm' m' 



dLuD'\ (lu)D''', {uj)D'\ (w) Stt' , , , , , 



where duj — sin fidpdad'j, to obtain 



-|-m 



£1 £2 i\f £1 £2 I 



^ \"'t^ra^^l2m2l ^ - l^va'^i2'nr,'^\ V M "^1 "^2 

.£1 ,mi ^2 ,*Ti2 m^m2 

K we now do the redefinitions 

, I (—1) 



and note that the first 3J symbol above is identically zero unless m[ — mj, we obtain finally (11 1. 



2. Useful identities 

We present here some useful identities related to the 3J symbols: 
• Isotropic limit 



mi 1712 J V2hTT " ' 



14 



• Parity and permutations 



Orthogonality 



mi 7712 ™ 



m mi 7712 



E E 

mi = — li m2=— /2 



E E (2^ + 1) 



(-1) 
(-1) 



h+h+i I h h I 

\ 7712 Wl Tn. 
/l+i2+i [ ^1 ^2 ^ 

-mi — m2 ~m 



h h h 
mi m2 ma 

'2 '3 
mi m2 m3 



E (-1)' 



m—~l 



ll I2 ^3 

mi m2 mg 

/i Z2 ^4 

mi m2 mg 

' / ; ^ 

m — m 



2/3 + 1 



The last expression is particularly useful in the derivation of (131. 



3. Some properties of the integral (12| 



The geometrical coefficients Igfg^ defined in (12 1 has many interesting properties which can be explored in order to 
speed up numerical computation of (11 1. First, we note that it is symmetric under permutation of £1 and £2 



E^t 



-ir 



ll £2 I 



+t2+l 



— E]^^2m,£im( 1) 



m+2{ei+e2+l) 



£2 h I 

—m m 



£2 ii I 

m —777 



Some of the other properties are a consequence of the integral defined in (24 1. We may note for example 

that, due to the symmetry of the X^m coefficient defined in (22 1, we will have 



leX^O, for any {(4, ^2) e N Ki + ^ = odd} . 

Furthermore, the Xgm coefficients restrict the m summation above to their values which obey m + £1 + £2 — even. 
K we further notice that ( 24 1 is proportional to the integral of a integral of the form Pe{cos 6) cos rnO d6, and that 
this integral is zero unless^ + m = even, we conclude that 

feX^O, for any {(^i, ^2, ^) e N ^i + 4 + ^ = odd} 
Besides, using the fact that the integral Pe {cos 9) cos m9 d9 is zero for any m < £, we find 



ri.e 



0, for any {{£i,£2,£) eN\£i < £,£2 < i}- 



We finally comment on the special case where I = 0, for which we have 
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However 

2f + l 



V2£ 



(2f + 1) (f + 



47r 



-Sfi' 



An {l' + m)\\ {l'-m)\\ 
Peix)Pe,ix)dx 



d{— cos -)?) , £' + m = even 



27r 

where in the derivation above we have made use of the Fourier series expansion of the Legendre polynomial. So we 
conclude that 



27rV2£i + 1 



(25) 



which is needed in the derivation of (13 1 and (18 1. 
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